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Abstract 

We observe signals for the decays ^(3770) — > XJ/ifj from data acquired with the CLEO detector 
operating at the CESR e + e~ collider with • v /s=3773 MeV. We measure the following branching 
fractions 6(^(3770) -»• XJ/tfj) and significances: (189 ± 20 ± 20) x 10~ 5 (11.6a) for X = vr+vr-, 
(80 ± 25 ± 16) x 10~ 5 (3.4a) for X = 7r°7r°, and (87 ± 33 ± 22) x 10~ 5 (3.5a) for X = n, where the 
errors are statistical and systematic, respectively. The radiative return process e + e~ — > jip(2S) 
populates the same event sample and is used to measure r ee [^(25)] = (2.54 ± 0.03 ± 0.11) keV. 
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The -0(3770) charmonium state decays most copiously into DD pairs, but other decays 
similar to those of 0(25) are predicted II S, B 0, 0, . The ip(2S) mass eigenstate is 
expected |8( have a dominant 2 3 Si angular momentum eigenstate with a small 1 3 D\ admix- 
ture, and vice versa for the ^(3770). Because more than half of ip{2S) decays contain a J /if) 
in the final state, a 2 3 S' 1 component enhances similar transitions for ■0(3770). Theoretical 
estimates (E IE IE 0] of the rate for transitions from the 1 3 Di eigenstate, based on a QCD 
multipole expansion, span a broad range. BES reported the first sighting of a 0(3770) non- 
DD decay 0, at ~ 3a significance, with £(0(3770) -> tt+tt- J/ip) = (0.34 ± 0.14 ± 0.09)%. 

In this Letter we describe a search for the XJ/ip final states, where X = tc + 7t~, 7r°7r°, 
r], and 7r°, in e + e~ collision data taken at a center-of-mass energy i/i=3.773 GeV. We use 
J /if) — > £ + £~, where £^ = or fj J ± . The data were acquired with the CLEO detector [lol ] 
operating at the Cornell Electron Storage Ring and correspond to an integrated lumi- 
nosity EE| of C = (280.7 ±2.8) pb _1 . The process e + e~ — > r yip(2S) dominates this event 
sample and is treated as background; it also yields a measurement of r ee [0(2S)]. 

The primary background for -0(3770) — ► XJ/ip is the tail of the 0(25) and radiative 
returns to it via initial state radiation (ISR) (i.e. e + e _ — > 7-0(25) — > -yXJ/ip); the total 
radiated energy peaks near 87 MeV but can take on a range of values. Similarly, there are 
radiative returns to that portion of the -0(3770) lineshape lying below ^/s which constitute 
part of the signal. The differential cross section for e + e~ — > ^R — > jXJ/if), where R = if>(2S) 
or -0(3770), can be expressed [3 EE EH in terms of XJ/ip mass-squared s' and the scaled 
radiated energy x = 1 — s ' / s as 

^ = W(s J x)xb(s , )xF x (s')xT ee xB x , (1) 

where W(s,ar) is the ISR 7-emission probability, b(s') is the relativistic Breit-Wigner for- 
mula, Fx(s') is the phase space factor |l2( appropriate for X, T ee is the e + e~ partial width 
of R (including vacuum polarization effects), and Bx = B{R — > XJ/ip) signifies an exclusive 
branching fraction. The ISR kernel is, at lowest order, 

W {8 i X>X ) = — (\n~l \ (l-x+^), (2) 
nx \ mi j \ 2 J 

in which xq > is a cutoff to prevent the divergence of / W dx, m e is the electron mass, and 
a is the fine structure constant. The Breit-Wigner function is 

in which Tr is the full width and Mr the nominal mass. The phase space factor is Fx{s') = 
(Px/Po) 2L+1 j i n which px is the momentum of X in the R center-of-mass frame, po is the 
value of px at ^fs' = Mr, and L is the relative orbital angular momentum between X and 
J /if). Eq. (1) has one enhancement near x = due to the 1/x factor in W(s,x), and, for 
s sufficiently larger than Mr, a much larger one near x = 1 — M R /s, corresponding to the 
peak of the Breit-Wigner resonance function. 

The cross section a(s) for e + e~ — > 7/? — > 7 XJ/ip can be both obtained from Eq. (1) and 
measured: 

= = F ee >< X , (4) 
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in which N is the number of events counted and e is the detection efficiency obtained from 
Monte Carlo (MC) simulation, and the integral 



I(s) = J W{s,x)b{s')F x {s')dx (5) 

is insensitive to the value of V. Hence a measurement of a(s) for R = ip(2S) can be combined 
with Bx measurements [lH to yield Y ee [%p{2S)}. 

We choose the E 1 cutoff to be 2 MeV (x = 1.06 x 1(T 3 ), small enough that events with 
x < Xq are experimentally indistinguishable from those with x = Xq. The expression 0,0, 
3 for W(s, x < Xq) includes terms accounting for soft and virtual photon emission (but not 
vacuum polarization, which is included in r ee ); we obtain I(s,x < x ) = 0.62 b(s) F x (s), a 
result reproduced by the Babayaga [l3[ [l + \i~ event generator. The integral I(s, x > Xq) can 
be performed numerically. For W(s, x > xq), we employ the full expression including higher 



order radiative corrections as given in Eq. (28) of Ref. 14]; it gives values ~19% smaller 
than Eq. (2) for J/ip radiative returns from y/s = 3.773 GeV. 

The EvtGen event generator [lj|> which includes final state radiation [2oj, and a GEANT- 
based j21| detector simulation are used to model the physics processes. The generator 
implements a relative 5*-wave (P-wave) configuration between the ttit (rj or 7r°) and the 
J/ijj. Radiative returns to ip(2S) and ?/>(3770) for x > x are generated with the polar 
angle distribution from Refs. jl6j], and account for ISR according to Eqs. (1-3). Separately, 
XJ/ip events are also generated without a photon to represent all x < xq events from 
ip{2S) and ^(3770), and are weighted with respect to the x > xq events according to the 
I(s,x < xo)/I(s,x > xo) ratios. 

Event selection implements the same requirements as in the CLEO ip(2S) — > XJ/ip 
analysis [3| except for the changes described here. No X-recoil mass cuts are imposed. 
To increase acceptance for J/ip — > £ + £~, lepton candidates at small polar angles (0.85 < 
|cos#£| < 0.93) are added. A number of measures are taken to reduce backgrounds. We 
demand m(£ + £~) = 3.05 — 3.14 GeV, and add to each lepton momentum vector any photon 
candidates located within a 100 mrad cone of the initial lepton direction. All 7r° candidates 
must satisfy 771(77) = 110-150 MeV. For X = ti + tt~ , neither pion candidate can be identified 
as an electron if m(7r + 7r~) < 450 MeV, which suppresses e + e~ — > £ + £~ r y, 7 — > e + e~ events in 
which the e + e~ pair from the photon conversion is mistaken for the 7r + 7r~. For X = 7t + tt~ 
and 7r°7r°, we require m(7nr) > 350 MeV. For the 7r°J/ip(—>- e + e~) andr/(— > 77)J r /'0(~ > e + e~) 
modes, background from Bhabha events is diminished by requiring cos 6 e + < 0.3. For X = 7r° 
or rj(— > 77), radiative transitions from ^(3770) or the ip(2S) tail to Xcj are suppressed by 
requiring the least energetic photon in the 7r° or rj candidate to satisfy > 280 MeV or 
£ 7 = 30-170 MeV. 

To extract the number of ij){2S) and -^(3770) events, we fit the distribution of event 
missing momentum, which can be interpreted as a measure of E 1 , 



s-Mj + m 2 x - 2Js(p x + m x 

k = -, — r 



2 (Jpi+Mj-pj cos. 



(6) 



in which Mj is the nominal J/ip mass, pj is the measured dilepton momentum, px is the 
measured X momentum, mx is the mass of X (the PDG value [22J for X = rj , 7r°, or the 
measured mass for X = tttt), and is the measured angle between the J/tp and the event 
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missing momentum three- vector (k). The small (~2 mrad) crossing angle of the incoming 
beams has been neglected. 

The phase space (l7[ factor for F n7T (s') is not as simple as the {px/pof dependence for 77 
and 7T° because the irn mass varies. The average momentum of the inr system increases by 
~11% from v / s~'=3.686 GeV to 3.773 GeV. As the nir and J/ip are in a relative S"-wave, we 
set F 7r7r (s'=(3.773 GeV) 2 ) = l.ll; for other s', F T7T is scaled linearly with x. The functional 
form of F nn is not crucial because da/dx is small over the central portion of the interval 
£ 7 = 0-87 MeV. 

The distribution in k for each exclusive mode is subjected to a maximum likelihood fit for 
three components with floating normalizations: a radiative return to ip(2S) shape obtained 
from MC simulation, a direct decay -0(3770) — ► XJ/ip signal shape from MC simulation 
(including radiative returns to the ■0(3770) tail), and a background component linear in k. 
Direct decays from the ^(3770) and the tail of the ip{2S) add incoherently 0, so that the 
ip(2S) background can be included without regard for interference. 

The distributions and fits are shown in Figs. ^El The fit results and quantities de- 
rived from them appear in Table The efficiencies include the MC correction factors from 
Ref. [lij . the visible cross sections use B(J/ip -> £+£-) from Ref. 0, and the T ee values 



use the B(ip(2S) — > XJ/ijj) results from Ref. jl8|. Statistical significances of the -0(3770) 
signals, obtained from the differences in log-likelihoods between fits with and without a 
signal component, are shown, indicating an unambiguous 7r + 7r~J/0 signal, and suggestive 
n n°J/ip and 77 J/-0 signals. The product of the measured cross section a(DD) 24] and lu- 
minosity C is used to give the number of produced -0(3770) decays as (1.80±0.03io!o2) x 10 6 - 
Sidebands around M(J/ip) in the dilepton mass distributions for events near the radiative 
return peak in k do not show evidence for additional background. Feed-across background 
(from radiative returns to ip(2S) but with non-signal il>(2S) decays) levels and uncertainties 
are determined using measured branching fractions 0] and MC simulation, as in Ref. [l8| . 

The results for 7-0(25') — > jit J/ip are all treated as upper limits due to substantial 
backgrounds from radiative Bhabha and muon pair events. Efficiency-corrected, background- 
subtracted rates for J/ip — > e + e~ and J/ip — ► n + fi~ are consistent with each other. Allowing 
second or third order polynomials in the background parameterizations has a negligible effect 
upon r ee . 

Statistical errors dominate for the 0(3770) results and systematic errors dominate for the 
ip(2S) results. Table HP summarizes the uncertainties that are uncorrelated for different X. 
The systematic errors on the fitted event yields include changes induced by variation of the 
range in k of the fit and the alternate use of a y 2 fit instead of maximum likelihood. The 
efficiency uncertainties are larger than in Ref. |l8[ because here the leptons are not restricted 
to | cos < 0.83, the 7r° — > 77 and J/0 — > £ + £~ mass cuts are tighter, and we account for 
imperfections in the assumed i/j(2S) boost direction. 

Relative uncertainties that are correlated for all X include those from B(J/ip — ► £ + £~) 
(0.94% statistical, 0.71% systematic), I(s) (i.e. radiative corrections) (2.0%), C (1.0%), and 
the normalization portion of Bx (3.0%). The MC sample used for this analysis has a mean 
and spread of very close to that of the data, within 0.05 MeV and 0.02 MeV, respectively, 
rendering negligible any remaining systematic effect upon I(s). 

A single value, Y ee [ip{2S)] = (2.54±0.03±0.11) keV, is obtained by combining TT + n~J/tp, 
7r 7r°J/?/>, and rjj/ip results, weighting each by the uncorrelated statistical and systematic 
errors. The relative 4.4% total uncertainty is dominated by the common 3.0% systematic 
normalization uncertainty in all CLEO ip(2S) branching fraction measurements ^j]. It 
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is 2.5 standard deviations higher than and of comparable precision to the PDG fit value, 
(2.12 ± 0.12) keV [i^; it is within two standard deviations of any of the results obtained 
from scanning the ip{2S) peak, the most precise of which is the preliminary BES [25| result, 
(2.25 ± 0.11 ± 0.02) keV. 

Figure |U shows the m(7r + 7r~) and £ + polar angle distributions for 7r + 7r~J/-0 events re- 
stricted to k = [—10, +10] MeV, background-subtracted with the sum of the [—57,-17] 
and [+13, +53] MeV sidebands scaled down by a factor of four. The MC histograms are 
normalized to the same areas as the data. In both plots, neither data nor MC is corrected 
for detection efficiency The data points represent events from both -0(3770) and ip{2S) in a 
ratio of ~2:1. The measured m(7r + 7r~) and \ cos6(£ + )\ distributions show consistency with 
the 0(2S')-like 5-wave MC predictions. 

The branching fraction for -0(3770) — > 7r + 7r~ J/tp is smaller than that reported by BES [9], 
but is consistent with it and more precise. While the widths for -0(3770) — > 7nrJ/ip are in 
the broad range predicted by the QCD multipole expansion models jHHEI 0; the 7171 mass 
distribution appears to be much stiffer than predicted for the large D-wave proportion fea- 
tured in these models. The branching fraction for ^(3770) — > ttttJ/i/j also relates to the 
interpretation of the X(3872): the small value does not strengthen the case for conventional 
charmonium |26[. The results combine to give £x£>('0(377O) — > XJ/ip) = (0.36 ± 0.06)%, 
which corresponds to a cross section of (23 ± 5) pb. The 90% C.L. upper limits are 
£>(-0(377O) — > rjj/ip, Ti°J/ip) < 0.15%, 0.028%; substantially more data would be required 
to quantitatively probe the cc purity of the ip{2S) and -0(3770) as proposed in Ref. |27| . 

In summary, we have observed the first statistically compelling signal for non-Dl) de- 
cays of the -0(3770), and with the same data sample have achieved improved precision on 
T ee m2S)}. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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TABLE I: Results for radiative return process e + e — > jip(2S), ip(2S) — > XJ/t/j and direct decay Y>(3770) — > XJ/tp. For each appears the 
fit yield AT, efficiency e, and cross section cr. In addition, for the radiative return process, the integral I(s) [followed by its value for x < xq], 
and the B{ip{2S) -► XJ/fp) x r ee values inferred from o appear along with the resulting r ee . The bottom five rows include the significance 
in standard deviations of the ^(3770) — ► XJ/tp signals and the ^(3770) branching fraction and partial width. Errors shown are statistical 
and systematic, respectively. 



X 


7T + 7T 


ttV 


V 


vr 


N(jtP(2S) -»• jXJ/iP) 


19469 ± 145 ± 195 


3616 ± 64 ± 72 


291 ± 19 ± 15 


< 37 


e(^{2S) - jXJ/i?) (%) 


56.23 ±0.07 ±0.90 


21.66 ±0.06 ±0.65 


7.89 ±0.17 ±0.28 


11.33 ±0.12 ±0.66 


a(^(2S) - 7 XJ/^) (pb) 


1036 ± 13 ± 23 


500 ± 10 ± 19 


111 ±8±8 


< 10 


J(s) (pb/keV) 


1215.4 [6.7] 


1215.4 [6.7] 


1251.9 [34.2] 


1215.2 [8.9] 


£(Y>(2S) X J/V) x T ee (eV) 


852 ±10 ±26 


411 ±8± 18 


88 ± 6 ± 7 


< 8 


r ee [^(2,S)] (eV) 


2541 ±32 ±113 


2488 ± 54 ± 138 


2716 ± 191 ±217 


< 6.2 x 10 3 


7V(Y>(3770) -► XJ/Y>) 


231 ± 24 ± 23 


39 ± 12 ±8 


22±8±6 


< 10 @90% C.L. 


Significance 


11.6(7 


3.4ct 


3.5a 


0(7 


e(^(3770) XJ/ip) (%) 


57.05 ±0.16 ±0.91 


22.86 ±0.13 ±0.69 


11.80 ±0.13 ±0.43 


16.02 ±0.15 ±0.93 


a(V>(3770) -► XJ/VO (pb) 


12.1 ±1.8 ±1.2 


5.1 ±2.0 ± 1.0 


5.5 ±2.1 ±1.4 


< 1.8 @90% C.L. 


£(^(3770) -► XJ/Y>) (10~ 5 ) 


189 ± 20 ± 20 


80 ± 25 ± 16 


87 ± 33 ± 22 


< 28 @90% C.L. 


T(V>(3770) XJ/VO (keV) 


45 ± 5 ± 7 


19±6±4 


21 ±8±6 


< 7 @90% C.L. 



TABLE II: Uncorrelated relative uncertainties in 
errors, see text. 



percent for the results in Table I; for correlated 



X 


7T + 7T 


A 


V 


7T° 


tp(2S) yield (stat) 


0.7 


1.8 


6.5 


19 


VK3770) yield (stat) 


10 


31 


36 




ij(2S) yield (sys) 


1.0 


2.0 


5.0 


10 


^(3770) yield (sys) 


10 


20 


25 




Efficiency 


1.6 


3.0 


3.6 


5.8 


i/)(2S) Feed-across 


0.1 


0.1 


3.0 


50 


B(ip(2S) -► XJ/il)) stat 


0.4 


0.9 


1.9 


7.7 


B{i){2S) -► XJ/ij)) sys 


1.3 


1.8 


1.6 


7.1 
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10 3 



i — i — | — i — i — i — r 
Jt + JT~J/lp 



T 1 1 T 

0* 



i — i — i — r 



2540805-012 
1 — i — i — i — - 



k (MeV) 




FIG. 1: Fit of the distribution in k for the final state tt + tt~ J/i/j, showing the data (open circles), 
overall fit (thin solid line), direct -0(3770) decay peak (thick solid line), radiative return to the 
ip(2S) (dotted line), and the background term (dashed line), on a logarithmic vertical scale (top), 
and on linear vertical scales focussed on the direct decay peak (bottom left) and radiative return 
peak (bottom right). 
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2540805-013 




k (MeV) 

FIG. 2: Fit of the distribution in k for the final state J /ip, showing the data (open circles), 
overall fit (thin solid line), direct V(3770) decay peak (thick solid line), radiative return to the 
ip(2S) (dotted line), and the background term (dashed line), on a logarithmic vertical scale (top), 
and on linear vertical scales focussed on the direct decay peak (bottom left) and radiative return 
peak (bottom right). 
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FIG. 3: Fit of the distribution in k for the final state rjJ/ip (top) and n°J/Tp (bottom), showing 
for each the data (open circles), overall fit (thin line), direct ■0(3770) decay peak (thick solid line), 
radiative return to the ip(2S) (dotted line), and the background term (dashed line). 
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FIG. 4: Distributions in tt + tt fi + l events of the tt + tt mass (left) and polar angle (right) of the 
positively charged lepton from data (open circles) and MC (solid line line). 
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